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^ , Abstract 



The experimental prospects for precise measurements of the leptonic decays B u — > 
rvjixv, B s — > /-i + /U _ , D — >■ fxu and D s — > \xvJtv are very promising. Double 

F-- ratios involving four of these decays can be defined in which the dependence on the 

values of the decay constants is essentially eliminated, thus enabling complementary 

^vq ■ measurements of the CKM matrix element V u b with a small theoretical error. We 

quantify the experimental error in a possible future measurement of \V u b\ using this 
approach, and show that it is competitive with the anticipated precision from the 
conventional approaches. Moreover, it is shown that such double ratios can be more 
effective than the individual leptonic decays as a probe of the parameter space of 
supersymmetric models. We emphasize that the double ratios have the advantage 

S^ ■ of using \V u b\ as an input parameter (for which there is experimental information), 

while the individual decays have an uncertainty from the decay constants (e.g. /b s ), 
and hence a reliance on theoretical techniques such as lattice QCD. 



1 Introduction 

Measurements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [1,2] are of 
great importance, being fundamental parameters of the Standard Model (SM). The matrix 
elements display a strong hierarchy, and the smaller off-diagonal terms can be expressed 
as higher powers of an expansion parameter A(~ 0.22) [3]. One of the smallest elements 
is V u b, being of order A 3 and with a magnitude ~ 4 x 10 -3 . Although the measurement 
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of \V u b\ is very challenging, in recent years its precision has been reduced to ~ 10% at 
the B factories. Two ways to measure \V u b\ have been utilized: inclusive semi-leptonic 
decays (B — y X u lv) and exclusive semi-leptonic decays (B — y iriv). These two ap- 
proaches provide independent measurements of \V u b\, which are currently not in total 
agreement, but have roughly the same precision. Both of these methods can be employed 
at high- luminosity B factories, and the experimental error in the determination of \V u b\ 
will decrease with increasing integrated luminosity and improved analysis techniques. The 
inclusive and exclusive determinations of \V u b\ have independent theoretical errors, and it 
is acknowledged that a concurrent reduction of this theoretical error (i.e. comparable in 
degree to that for the experimental error) is more challenging (e.g. see [4]). Therefore, as 
emphasized in [5-8] it is important to consider alternative methods of measuring \V u b\ for 
which the theoretical error is known to be very small, even if the experimental error for 
such a method is currently large. 

One such approach is the branching ratio (BR) of the purely leptonic decay B u — y 
tv. The unprecedented data samples provided by the B factories have enabled the first 
measurements of B u — y tv [9-12] despite its relatively small BR and challenging signature. 
The current world average of BK(B U —y tv) has an error of 20%, although proposed 
high-luminosity flavour factories could reach a precision of 5% [13-21]. A measurement 
of |Kif,|/s can be obtained from this decay, where /# is the decay constant for the B u 
meson, which has been calculated by Lattice Quantum Chromodynamics (LQCD). An 
average of two unquenched calculations [22,23] of /g gives /# = 192.8 ± 9.9 MeV [24], 
i.e. a precision of ~ 5%. Hence the decay B u — y tv (and its companion decay B u — y fiv, 
which is yet to be observed), could provide complementary (and possibly competitive) 
measurements of \V u b\ at high-luminosity flavour factories. 

However, despite the impressive improvement in the precision of calculations of fs 
by LQCD, it is important to extract \V u b\ with minimal theoretical error (i.e. with 
minimal dependence on parameters like /#, which are determined by the strong in- 
teraction). It was pointed out long ago [5] that a double ratio of the leptonic decays 
B u — y tv,B s —y /u + fi~,D —y \xv and D s —y \wjrv has essentially no dependence on the 
absolute values of the decay constants. The decay constants in this double ratio appear 
as (/b s //b)/(/d s //d), and such a quantity is numerically very close to 1, with a very 
small error. Hence such a double ratio can be used to obtain a measurement of |V^&| with 
a small theoretical error [7]. At the time of writing of [5], only one of the four decays had 
been measured (D s —y \ivJtv). At present, three of the four decays have been observed, 
and the prospect of precise measurements of all four decays in this double ratio motivates 
us to quantify the expected precision for a measurement of \V u b\ using this method. 

The decays D ->■ pi/ [25], D s -> fiv [26-28] and D s -> tv [26,29-31] have been 
measured with ~ 10% precision at CLEO-c and the B factories, and the ongoing BES-III 
experiment will reduce the error in all these decays to maybe 3% [32,33]. The current 
measurements of B u — y tv, D —y fiv, D s — y fiv and D s —y tv are in agreement with the 
SM prediction, within the errors of the decay constants. An upper limit exists for the 
decay B s —y n + fi~ (< 4.3 x 10^ 8 at 95 C.L. [34]) and this is the only decay in the double 
ratio which has not been observed. At the CERN LHC, the LHCb experiment is expected 



to achieve sensitivity to the SM prediction of BR(5 S — > fi + fi~) ~ 3 x 10~ 9 . Therefore the 
first measurement of the double ratio will be possible by combining the measurements of 
these four leptonic decays from three distinct experimental facilities. 

The above purely leptonic decays are sensitive to New Physics particles which arise 
in supersymmetric (SUSY) models, such as charged Higgs bosons (i? ± ), neutral Higgs 
bosons (h°, H° and A ), and charginos (x ± )- The decays B u — > tv, D — > fiv, D s — > \iv and 
D s — > tv are sensitive to H ± at the tree level [35]. The largest effect is for B u — > tv and the 
current measurement of this BR, although in agreement with the SM prediction, does not 
preclude a large contribution from H^ [36] 3 . The effect of H ± on D s — > fiv and D s — > tv 
is much smaller [35,37], but can be of the order of the current experimental precision 
[38-40]. The contribution of H ± to D — > fiv is negligible. The decay B s — > fi + fi~, which 
proceeds via higher-order diagrams in perturbation theory, can be greatly enhanced by 
contributions from neutral Higgs bosons and chargino-stop loops. The contributions from 
SUSY grow with tan/3, which can lead to an enhancement of BK(B S — > fi + fi~) by an 
order of magnitude with respect to the SM expectation [41,42]. 

Since the decays involved in the leptonic double ratios can receive substantial con- 
tributions from Higgs bosons and supersymmetric particles, another possible application 
of the double ratio is to constrain SUSY parameters. The absence of an explicit depen- 
dence on the decay constants makes the double ratio a particularly attractive quantity 
with which to probe SUSY. In this case, V u b is an input parameter and is the only major 
source of uncertainty. Information on V u b is available from both direct measurements 
(the aforementioned inclusive and exclusive semileptonic decays, the latter having some 
dependence on LQCD) and from global SM fits. In contrast, the calculation of decay con- 
stants is reliant on LQCD. Consequently, due to its reduced theoretical uncertainties the 
double ratio constitutes a cleaner observable than the individual decays, and can provide 
competitive constraints on SUSY parameters. 

The paper is organized as follows: in section 2 the conventional ways to measure \V U b\ 
are briefly summarized and the double ratio of leptonic decays is introduced. In section 
3, the decay rates of the leptonic decays are presented and the experimental prospects for 
their measurements are discussed. The current upper limit on the double ratio is derived, 
as well as the anticipated error in a measurement of \V u b\ at future experimental facilities 
using the double ratio. Section 4 contains a numerical study of the double ratio as a probe 
of the SUSY parameter space in two distinct models. Conclusions are given in section 5. 

2 Measurements of \V u b\ 

Currently there are two methods which are used to measure \V u b\- Both methods use a 
semi-leptonic decay channel of the B/B u meson for which the decay rate depends on | V u b\ 2 . 
One method ( "inclusive approach" ) is to measure the branching ratio of the inclusive decay 
B — > X u £v while the other method ("exclusive approach") is to measure the branching 
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ratio of the exclusive decays B — > ir£v and B — V p£v. The value of \V u b\ is then extracted 
from these measurements, with an additional error from the theoretical calculation of 
the decay rate, which is independent for the two decays. Currently there is a small 
disagreement between the values of \V u b\ extracted from these distinct approaches. 

2.1 Inclusive determination of |T4&| 

The main difficulty with the inclusive determination of \V u b\ is the fact that the decay 
B — > X c £v has a branching ratio about fifty times larger than that of B — > X u £v (because 
\V c b\ 3> |T4&|)- In order to reduce this background and isolate the signal from B — > X u £v, 
kinematical cuts are applied. One approach is to select charged leptons £ with greatest 
energy (the "endpoint region"). For an energy Eg_ > 2.4 GeV the background from 
B —?• X c £v is very small, but the theoretical prediction for the decay rate of B — > X u £v in 
this endpoint region has sizeable errors, which are summarized in [4] . The most significant 
error is from the non-perturbative "shape function" [43]. Fortunately, the leading shape 
function is universal in B jB u decays. Therefore the dependence of the decay rate of 
B — > X u £v (in the endpoint region) on the leading shape function can be eliminated 
by using data from other B decays (such as B — >• X s j). However, "sub-leading shape 
functions" break this universality of the shape function, and so the process of eliminating 
its effect on B — > X u £v brings in additional errors. In [8] it is suggested that this error 
from sub-leading shape functions might limit the ultimate precision in such measurements 
of \V u b\ to around 15%. Alternatively, the dependence of the decay rate of B — > X u £v on 
the shape function can be reduced by applying a less energetic cut on £ (e.g. E( > 1.9 
GeV), but the background from B — > X c £v then becomes much larger. 

There are several measurements of | V^t, | from the inclusive method, and different kine- 
matical cuts are applied. In order to reduce the background from B — > X c £i>, most 
measurements select only a small fraction (partial rate) of the total inclusive rate of 
B —?• X u £v. The value of \V u b\ is then extracted in the context of several theoretical 
prescriptions, which treat the various theoretical errors (including the dependence on the 
shape function) in distinct ways. 

A recent measurement by the BELLE Collaboration [44] applies a much lower cut 
on the energy of £ (Eg > 1.0 GeV), which keeps ~ 90% of the total inclusive rate of 
B — > X u £u. Such an approach substantially reduces the theoretical error (especially 
that from the shape function) in the extraction of \V u b\, and it is suggested in [16] that 
this method is the most promising one for a precise measurement of \V u b\ from inclusive 
B — > X u £v at a high-luminosity B factory. 

The current measurement of \V u b\ from inclusive decays uses eight independent mea- 
surements of B — \ X u £v, and is an average over the results of the four theoretical pre- 
scriptions [45]: 

|K&| =4.37 ±0.39 x 10~ 3 . (1) 

In Eq.(l) the theoretical error and experimental error have been combined, and the ex- 
perimental error (±0.16) is smaller than the theoretical error. 



2.2 Exclusive determination of \V u b\ 

The decays B — > nlv and B — > plv have both been measured, but the former is much 
more favourable for extracting \V u b\, both experimentally and theoretically. There are two 
search strategies for B — > nil/. One technique tags the other B/B u in the event, while 
the other technique does not. Currently these techniques have comparable precisions, but 
it is suggested that the tagged method (e.g. as used by BELLE in [46]) will afford the 
greatest precision at a high-luminosity B factory. 

Extracting the value of \V u b\ requires theoretical input distinct from that used in the 
inclusive approach. The decay B — > iriu only depends on one hadronic form factor f+(q 2 ) 
(where q is the momentum of £) which can be calculated in LQCD for q 2 > 16 GeV 2 
and by Soft Collinear Effective Theory (SCET) for q 2 < 16 GeV 2 . The world average 
measurement of \V u b\ from exclusive decays is [24]: 

\V ub \ =3.42 ±0.37 x 10~ 3 . (2) 

A recent BABAR measurement alone [47] gives \V u b\ = 2.95 ± 0.31 x 10 -3 . Therefore 
both the inclusive and exclusive measurements of \V u b\ have comparable experimental and 
theoretical errors, but there is a disagreement of the central values. However, the average 
of the inclusive and exclusive measurements is \V u b\ = (3.92 ± 0.09 ± 0.45) x 10~ 3 [48], 
which agrees quite well with the values of [V^l extracted from two independent global 
SM fits, the latter being |K&| = 3.53±0.15 x 10~ 3 [48] and |Kb| = 3.52±0.11 x 10~ 3 [49]. 

2.3 The double ratio (/b//b.)/(/d//d.) 

Although the experimental prospects for both the inclusive and exclusive approaches 
to measure |V^,| are good, the question of the magnitude of the theoretical error will 
remain. A complementary technique for which the theoretical error is known to be small 
is particularly attractive. In this regard, the use of double ratios of leptonic decays to 
measure \V u b\ was proposed in [5,7]. One such double ratio is defined by: 
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T(B S -> fi+pi-) T(D -> H |^ fc | 2 a 2 f D /f Da " 

The quantity (/b//bJ/(/d//d s ) deviates from unity by small corrections of the form 
m s /mb and m s /m c , and importantly such a double ratio does not have a dependence on 
the absolute values of the decay constants. A calculation in [5] gives (/s//sJ/(/d//d s ) = 
0.967, and subsequent works [50] also give values very close to 1, with a very small error. 
In our numerical analysis we will take (/b//b s )/(/£»//d s ) = 1- Therefore Eq. (3) can 
be used to measure \V u b\ with essentially no dependence on non-perturbative techniques 
such as LQCD. This amounts to using the decay B u — > rv to measure |T4&|, with the 
dependence on /# (i.e. the dependence on non-perturbative techniques) being removed 
by forming the double ratio of decay constants. In Eq. (3), the product of CKM elements 
I Vt s Vtb\ 2 is known with high precision if one assumes that the CKM matrix is unitary. The 
left-hand side of Eq. (3) is taken from experiment, and thus a measurement of \V u b\ can 



be extracted when all four decays have been measured. The drawback of this approach 
is that accurate measurements of all four decays are required. At the time of writing 
of [5] only one decay had been measured. At present, three of the decays in Eq. (3) 
have been measured, while an upper limit exists for B s — > fi + fi~. Evidently, the double 
ratio cannot yet provide a measurement of \V u b\ because one of the decays has not been 
measured. However, the experimental prospects for reasonably precise measurements of 
all four decays are quite good, and eventually a measurement of \V u b\ using Eq. (3) will 
be possible. 

Other double ratios can be defined, such as replacing D s — > fiu in Eq. (3) by D s — > tv, 
both of these decays having comparable experimental precision at present (and similar 
precision is expected for both decays at the BES-III experiment [33]). Moreover, the 
decay B u — > fiis is expected to be observed at a high-luminosity B factory, which would 
provide an independent measurement of \V u b\fB, and enable a measurement of a double 
ratio in which B u — > tv is replaced by B u — > fiv. We note that a SM-like rate for the 
decay D — > tv will be measured at the BES-III experiment, but its precision will always 
be inferior to that for D — >• fiv [51], and so it is preferable to use this latter decay in the 
definition of double ratios. 

In the context of models beyond the SM, the double ratio is sensitive to New Physics 
particles. Assuming that the partial widths for each decay are multiplied by scale factors, 
the right-hand side of the double ratio in Eq. (3) is modified by a ratio of scale factors: 

T -^±- = R- 1 . (4) 

r B a 

In the next section we will give the explicit form of each of the scale factors appearing in 
R, in the context of the Minimal Supersymmetric Standard Model (MSSM). 

3 The leptonic decays in the double ratio 

In following subsections we present the theoretical expressions for the decay rates of 
the leptonic decays in the double ratio, and discuss the experimental prospects for their 
measurements. The current upper bound on the double ratio is derived by combining 
statistically the existing measurements of the individual decays. Finally, the attainable 
experimental precision for \V U b\ from a measurement of the double ratio at the next gen- 
eration of experiments is estimated. 

3.1 The decay B u — > tv 

In the SM this decay is mediated by W ± and is helicity suppressed, whereas there is no 
such suppression for the scalar H ± exchange, whose contribution is proportional to the b 
quark and r lepton Yukawa couplings. In the limit of high tan /3 such Yukawa couplings 
are enhanced, and the contribution from H ± can be comparable in magnitude to that of 



W ± [35]. The leading order SM prediction for this decay is: 
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while the New Physics contribution from i7 ± is expressed through the ratio [35, 52, 53] 
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Here rriH+ is the mass of the charged Higgs boson, m# is the mass of the B u meson, 
T B is the total decay width of the B u meson, tan/3 = v 2 /vi where V\ and v 2 are the 
vacuum expectation values of the two scalar doublets, and Gp is the Fermi constant. The 
leading SUSY-QCD corrections are included in this expression through 6q. The Yukawa 
couplings of the MSSM take the form of a 2HDM (Type II) at tree level, but at higher 
orders the structure becomes of the type 2HDM (Type III) in which eo is a function of 
SUSY parameters [54-57] and |eo| can reach values of order 0.01. Using the average f B = 
192.8 ± 9.9 MeV [24] of two unquenched lattice QCD calculations of f B [22,23], and the 
average of the inclusive and exclusive determinations of \V u b\ = (3.92 ±0.09 ±0.45) x 10~ 3 
[48], we evaluate the SM branching ratio with Superlso v2.8 [58,59]: 
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There are four independent measurements of B u — > tv [9-12] and the current world 
average is [60]: 

BR(5 U -¥ ru T ) cxp = (1.63 ± 0.39) x 10" 4 . (8) 

The SM prediction can be compared to the experimental average by forming the ratio: 
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(9) 



Significant improvements in the precision of the measurement of B u — > tv will require 
a high- luminosity B factory. The two measurements by BABAR [11, 12] have used a 
large amount (~ 70%) of the available data taken at T(4S). The two measurements 
by BELLE [9, 10] are with 414 fb _1 and 605 fbr 1 , and so these could be significantly 
updated with the total integrated luminosity of 1000 fb" 1 . A high-luminosity B factory 
with around 50 ab _1 could measure B u — > tv to a precision of around 6% (e.g. see [19]). 
Moreover, with 50 ab _1 the decay B u — > fiu (for which there is currently an upper limit) 
could be measured with about the same precision as B u — > tv. Hence a precision of 
around 3% for |T4&|/b could achieved from each decay at a high- luminosity B factory. 
A combination of the measurements of |VJi&|/£ from B u — > tv and B u — > fiv would 
(presumably) further reduce the uncertainty. 



3.2 The decays D s — > tv, D s — > \iv and D — > \w 

In analogy to the case for B u — > tv, singly charged Higgs bosons would also contribute to 
the decays D s — > tv and D — > tv at tree level [35] . The effect is negligible for D — > tv, but 
it can be of the order of the current experimental precision for D s — > fiv and D s — > tv. 
The partial width is given by (where t = e, \i or r): 



T(D a -+ lv t ) = ^flm]M Ds (l - ^-) \V cs \\ Da , (10) 
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where in the MSSM one has [35,37,39,40,61]: 
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Here m c and m s are the masses of the charm and strange quarks respectively, Md s is the 
mass of the D s meson, V cs is a CKM matrix element, and me is the lepton mass. We note 
that eo in Eq. (11) is not the same as eo in Eq. (6), because they are functions of different 
SUSY parameters. The term m s tan 2 (3, which originates from the strange quark Yukawa 
coupling, can give rise to a non-negligible suppression of ro s for large values of tan fl. 
Note that the magnitude of the H ± contribution depends on the ratio of quark masses 
m s /(m c + m s ), and an analogous uncertainty is not present for the H ± contribution to 
the decay B u — > tv. However, there is now a very precise calculation (error ~ 1%) of 
m s /m c from unquenched LQCD [62]. 

There are various unquenched lattice calculations of fn s and the current situation 
is summarized in [63]. The value with the smallest quoted error is from the HPQCD 
Collaboration, f Ds = 241 ± 3 MeV [64] (with a provisional update of f Ds = 247 ± 2 MeV 
[65]). The MILC collaboration obtains fo s = 249 ±11 MeV [66] (with a provisional 
update of fo 3 = 260 ± 10 MeV [67]). A partially quenched calculation by the ETMC 
Collaboration gives f Da = 244 ± 8 MeV [68]. 

On the experimental side, D s — > tv has been measured at CLEO-c for three decay 
modes of r (evv, nv, pv) [26,29,30]. These are absolute branching ratio measurements. 
Moreover, a measurement of D s — > tv has been performed at BABAR [31] for one decay 
mode of r (evv). Unlike the CLEO-c measurements above, this measurement is normal- 
ized to the branching ratio of the decay D s — > K ± K°. The decay D s — > fiv has been 
measured at CLEO-c [26], BELLE [27] and BABAR [28]. Both the CLEO-c and BELLE 
measurements are of the absolute branching ratio, while the BABAR measurement is 
normalized to the decay D s — > 0°7r ± (note that this differs from the decay D s — > K^K 
which is used in the BABAR measurement of D s — > tv). Due to the sizeable uncertainty 
in the branching ratio of D s — > 0°7r ± , the BABAR measurement [28] is not included in 
the averages given in [63]. The average of the six measurements [26,27,29-31] results 
in a world average of fu g = 257.5 ±6.1 MeV (derived in [63] and is dominated by the 
average of the four CLEO-c measurements). The world average given in [60] is lower 
(fo s = 254.6 ± 5.9 MeV) because it includes the BABAR measurement of fo a in [28] 
(although with a reinterpreted value, as explained in [69]). 
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The experimental results for the branching ratios are [40,60]: 
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(5.38 ±0.32) x 10" 



BR(L> S -»■ /^) cxp = (5.81 ± 0.43) x 10" 3 
BR(D -> /i^) cxp = (3.82 ± 0.33) x 10~ 4 



(12) 

(13) 
(14) 



More precise measurements will be possible at BES-III. A precision of a few percent 
(2% — > 4%) is expected for all three decays D s — > tu t , D s — > fiu^, and D — > /iv^, before a 
high-luminosity B factory starts to operate. In principle, data from the current B factories 
(and at high-luminosity upgrades) could also be used to provide further measurements 
of D s — > fiu and D s — > tv. The recent measurement of D s — > tv (with r — > evv) by 
BABAR [31] suggests that a similar measurement could be attempted at BELLE. Other 
decay modes of r (e.g. tt ± v and p ± v, as done at CLEO-c) might also provide additional 
measurements. We note that BES-III has the capability to measure the decay channel 
D s — > K^K (which is used for normalizing in [31]) more precisely, which will reduce the 
error in the extracted branching ratio of D s — > rv T from techniques like that used in [31]. 
The SM predictions for the BRs, obtained using Superlso v2.8, are as follows: 



BR(D S ->■ tv t ) = (4.82 ± 0.14) x 10~ 2 

BR(D S ->■ \xv^ = (4.98 ± 0.14) x 10~ 3 

BR(D -> /iz/ M ) = (3.89 ± 0.16) x 10~ 4 



(15) 
(16) 
(17) 



in which f Da = 241 ± 3 MeV is used. 

3.3 The decay B s -¥ fi + n~ 

The decay B s — > fi + fi~ is the only decay in the double ratio which has not been measured. 
In the SM, the main theoretical uncertainty is from /# s and V ts , but the latter can be 
taken as a well-measured parameter by assuming CKM unitarity. It has been emphasized 
in many works [41,42,70-72] that the decay B s —$■ ji + ji~ is very sensitive to the presence 
of SUSY particles. At high tan/3, the MSSM contribution to this process is dominated by 
the exchange of neutral Higgs bosons. We therefore expect indirect constraints on m#+ 
and tan/3 from the MSSM mass relations. The BR(£? S — > fi + fi~) can be expressed as [73] 
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where the coefficients Cs, Cp, and Ca parametrize different contributions. Within the 
SM, Cs and Cp are small, whereas the main contribution entering through Ca is helicity 
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suppressed. In the MSSM, both C$ and Cp can receive large contributions from scalar 
exchange. The B s decay constant fg a = 238.8 ±9.5 MeV [24] constitutes the main source 
of uncertainty in this expression. 
We define 

BR(I? S — > ^ + /i")MSSM /lri N 

The SM prediction, obtained with Superlso v2.8 is 

BR(B S ->■ ^ + /i") S M = (3.21 ± 0.29) x 10~ 9 , (20) 

while the current experimental limit, derived by the CDF collaboration, is [34,60]: 

BR(B S -> fx + fi~) < 4.3 x 10" 8 (21) 

at 95% C.L. The experimental limit is thus still an order of magnitude away from the 
SM prediction, allowing for a substantial SUSY contribution. Including theoretical un- 
certainties, we compare the MSSM prediction to the upper limit at 95% C.L. 

BR(B S -► //+/0 < 4.7 x 10" 8 . (22) 

Three experiments at the LHC (LHCb, CMS and ATLAS) will have sensitivity to 
BR(i? s — y fi + fi~) which is superior to that at the Tevatron. At LHCb, a 3a signal would 
be established for a SM-like rate for BR(5 S — > /i + /i~) after one year at design luminosity 
(2 fb" 1 ) [74]. Slightly inferior sensitivity is expected at CMS and ATLAS [75], and a 3cr 
signal would require three years at design luminosity (30 fb _1 ). The dominant systematic 
error will be from normalizing to specific decays of B u and B, with the two main sources 
of error arising from: i) the uncertainty in the relative production rate of B s and B u /B; 
ii) the error in the measurements of the branching ratios of B u /B which are used as 
normalization channels. In [76] this error is estimated to be 13%, thus making it difficult 
to establish (at a high confidence level) enhancements of up to a factor of 3 compared 
to the SM expectation. An error of ~ 13% (not including statistical error) seems to be 
realistic for a measurement of a SM-like BR(i? s — > fi + fi~). 

3.4 Experimental limit on the double ratio 

The double ratio in Eq. (3) can be rewritten as 

R = — , (23) 

??SM 

where 

= / BR(B,->/iV) \ / ( BR(D.^tv) \ 

1 ~ V BR(B U -+tv) ) I V BR(£> -)• fiu) ) ' l ] 

As explained in section 2.3, this ratio does not suffer from uncertainties from the decay 
constants, and the main theoretical uncertainty in the evaluation of 7/sm (= (2.27±0.54) x 
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10~ 7 ) comes from V u b- To determine the experimental limit on the ratio (24), we combine 
the limits on the individual branching fractions in Eqs. (8), (12), (14) and (21). We 
use a Gaussian distribution for the constraints in Eq. (8), (12) and (14), and a "flat" 
distribution for the upper limit in Eq. (21). Combining statistically these experimental 
limits, we obtain, at 95% CL. 

7] < 2.09 x 10" 6 , (25) 

which yields the upper limit for R at 95% C.L. 

R < 10.0 , (26) 

in which the uncertainty from V ub is taken into account. In our numerical analysis we use 
Eq. (26) to constrain the supersymmetric parameter space in two scenarios of the MSSM. 

3.5 Attainable precision for \V u b\ using the double ratio 

It is possible to estimate the future error in a measurement of \V u b\ using the double ratio, 
assuming SM-like measurements for all four decays. The uncertainty in \V u b\ will come 
from the error in the measurements of the four leptonic decays (i.e. the left-hand side of 
Eq. (3)). 

We first use the current experimental errors for each measured branching ratio, and 
consider a SM-like value for BR(I? S — > fi + n~) with 30% relative experimental error. 
Calculating \V u b\ with a random Monte-Carlo generator, we find 

|K&| = (3.58 ± 0.71) x 10" 3 , (27) 

which corresponds to about 20% relative error. Although the error in Eq. (27) is larger 
than that for the inclusive and exclusive measurements in Eq. (1) and Eq. (2), the dou- 
ble ratio would provide an additional direct measurement of \V u b\ with reasonably good 
precision. Moreover, such a measurement of \V u b\ would be complementary, because the 
theoretical error in Eq. (27) is much smaller than that in Eq. (1) and Eq. (2). 

The future measurements of the branching ratios involved in the double ratio will have 
reduced errors, such as 2.5% for D and D s decays, and (in the optimistic case) 5% for 
BR(5 —} tu), and 13% for BK(B S — ¥ /z + /i~). With this set of errors, the uncertainty in 
the determination of \V u b\ would be reduced: 

\V ub \ = (3.58 ± 0.27) x 10" 3 , (28) 

corresponding to 7.5% relative error. Hence the double ratio has the potential to provide 
a measurement of \V u b\ which has a precision comparable (or even superior) to that for 
the current measurements of \V u b\ in Eq. (1) and Eq. (2). Moreover, the analogous double 
ratio with the decay B u — » \iv would lead to an independent determination of \V ub \, 
with an error which is comparable to that obtained from the double ratio with BK(B — > 
tv). Combining the measurements of these two double ratios could further improve the 
precision of \V u b\. 
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We note that the inclusive and exclusive determinations of \V u b\ will also have reduced 
experimental errors at a high- luminosity B factory [16], although the extent of the im- 
provement of the theoretical error is not clear at present. The inclusive determination of 
\V u b\ used in [44] is a very promising approach because it has a relatively small theoretical 
error, and so such a method is likely to lead to the most precise measurement of \V u b\ 
at a high-luminosity B factory. However, even if the total error from the conventional 
determinations of \V u b\ is smaller than that from the double ratios, the complementary 
measurements of \V u b\ from the latter will be important additional information, especially 
if the inclusive and exclusive determinations continue to give significantly different central 
values (which is the case at present). 

4 Constraints on the supersymmetric parameter space 

In this section we study the dependence of the double ratio R on the parameter space 
of supersymmetric models. We consider two scenarios in the MSSM: i) the constrained 
MSSM (CMSSM) and ii) the Non-Universal Higgs Mass model (NUHM), with minimal 
flavour violation (MFV). The CMSSM is characterized by a set of universal parameters 
at the GUT scale {mo, mi/2, ^4o 5 tan/3}, as well as the sign of the \i parameter, while in 
the NUHM the universality of the GUT scale mass parameters is relaxed for the Higgs 
sector leading to two additional parameters, \i and tjia- This additional freedom implies 
in particular that the charged Higgs mass can be considered as a free parameter. 

To investigate the parameter spaces of the CMSSM and NUHM we generate 300,000 
random points scanning over the ranges m G [50, 2000] GeV, myz G [50, 2000] GeV, 
A G [-2000, 2000] GeV, tan/3 G [1, 60] with positive // in CMSSM and // G [-1000, 2000] 
GeV and tua £ [5, 2000] GeV in NUHM. For each point we calculate the spectrum of 
SUSY particle masses and couplings using SOFTSUSY 3.1 [77] and we compute the 
branching fractions BR(B U — > tu t ), BR(D — ¥ fiv^), BR(D S —} ru T ), BR(D S — > \xv^) and 
BR(B S —?• /U + /i~) using Superlso v2.8 [58,59]. The obtained values are then used to 
calculate the double ratio. 

In Fig. 1, the results are shown for CMSSM in two separate plots. The scan over 
the four dimensional parameter space is projected into the plane (tan/3, min), where 
different colours correspond to different intervals for R. The points in black are excluded 
by R. In the left plot, the excluded points are displayed in the background, while they 
are shown in the foreground in the right plot. The excluded region in the left plot is 
therefore independent of the other SUSY parameters. In a large part of the parameter 
space, the double ratio is SM-like. In these regions, \V u b\ can be determined with almost 
no additional deviation due to SUSY. On the other hand, the area with tan/3 > 55 and 
mi/2 < 1000 GeV is excluded with no dependence on the lattice inputs. 

In a similar way, in Fig. 2 the results are shown for NUHM in the plane (tan^m^). 
Again, the points in black are excluded by R, in the background in the left plot and in 
the foreground in the right plot. As in CMSSM, in a large part of the parameter space, 
the double ratio is SM-like, while the region m^/tan/3 < 8 GeV is excluded with no 
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Figure 1: Variation of the double ratio R in the CMSSM plane (tan ft, mi/2). The zones 
in green and red delimit 1% and 10% deviation from the SM value respectively. In the 
yellow zone, R can be a factor of 2 away from the SM and in the grey zone by a factor of 
10. The black points are excluded at 95% C.L. 
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Figure 2: Variation of the double ratio R in the NUHM plane (tan/3, m^). The zones 
in green and red delimit 1% and 10% deviation from the SM value respectively. In the 
yellow zone R can be a factor of 2 away from the SM and in grey zone by a factor 10. 
The black points are excluded at 95% C.L. 



dependence on the lattice inputs. 

In order to compare the constraints on the SUSY parameter space from the double 
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Figure 3: In the left panel, CMSSM parameter plane (mi/2,m ) with tan/3 = 35, A = 
and /i > 0. In the right panel, NUHM parameter plane (/x, ttia) with tan f3 = 35, A = 
and rrtQ = mi/2 = 1000 GeV. The points in red are excluded at 95% CL. by BK(B S — > 
fi + [i~) and in grey by R. 



ratio with those from BK(B S — > /i + ii~) alone, in Fig. 3 we show examples of parameter 
space in the CMSSM and in the NUHM which are excluded from both observables. It 
can be seen that the double ratio, being a combination of four different constraints, 
can be more constraining than the branching ratio of B s — > /i + ii~ taken individually. 
Importantly, contrary to B s — > /i + /i~, the double ratio does not depend on lattice inputs, 
but instead depends on |V^,|, whose magnitude is already constrained from various distinct 
experimental methods (i.e. direct measurements and global-SM fits). Hence the double 
ratio is an important alternative to BR(5 S — > fi + fi~) as a probe of the SUSY parameter 
space. 

The conclusion that the double ratio can be more powerful at constraining SUSY 
parameters than BR(5 S — > fi + fi~) is irrespective of the freedom in choosing the input 
parameters (which are \V u b\ for the double ratio and fs s for BR(£? S — > /i + /i - )). This is 
shown in Fig. 4 where both the least favourable (with high \V u b\ and high fg s ) and most 
favourable (with low \V u b\ and low fg s ) cases for the double ratio are considered. In both 
cases, the double ratio excludes a greater region of the SUSY parameter space. 



5 Conclusions 

Ongoing and forthcoming experiments will significantly improve the precision of the mea- 
surements of the leptonic decays B u — > tv (at a high-luminosity B factory), D —¥ fxu 
(at BES-III) and D s — > [ivjrv (at a high-luminosity B factory and BES-III). A first 
measurement of the decay B u — > \iv will be possible in the early stages of operation of 
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Figure 4: In the left panel, CMSSM parameter plane (mi/2,m ) with tan/3 = 35, A = 
and \i > for \V ub \ = 3.48 x 10" 3 and f Bs = 238.8 MeV. In the right panel, \V ub \ = 4.5 x 
10- 3 and f Bs = 250 MeV. The points in red are excluded at 95% C.L. by BR(B S ->• ii + ^~) 
and in grey by R. 



a high-luminosity B factory if its branching ratio is comparable to the prediction in the 
Standard Model. The ongoing experiments LHCb, CMS and ATLAS will provide the first 
measurements of the decay B s — > fi + fi~ if its branching ratio is comparable to (or greater 
than) the prediction in the Standard Model. The branching ratio of each of these leptonic 
decays depends on the magnitude of the relevant decay constant, which can be calculated 
with techniques such as lattice QCD. Double ratios involving four of the decays can be 
defined in which this dependence on the decay constants is essentially eliminated, thus 
giving rise to an alternative observable with substantially reduced theoretical error. 

A measurement of the double ratio involving B u — > tu, D — > fiu, D s — > [ivjrv and 
B s —> H + fi~ would enable a complementary measurement of the CKM matrix element V ub 
[5,7] with a theoretical error which is much smaller than that present in the conventional 
approaches to measure \V u b\ (which utilize inclusive and exclusive semi- leptonic decays 
of B mesons). We quantified the experimental error in a possible future measurement 
of \V u b\ using the above double ratio, and showed that it can be competitive with the 
anticipated precision from the conventional approaches. Such an additional measurement 
of \V u h\ would be particularly beneficial if the current disagreement between the inclusive 
and exclusive determinations of \V u b\ persists into the era of a high- luminosity B factory. 

In the context of supersymmetric models the above leptonic decays are also mediated 
by New Physics particles. We showed that such double ratios can be more effective 
than the individual leptonic decays as a probe of the parameter space of supersymmetric 
models. We emphasized that the double ratios have the advantage of using \V u b\ as 
an input parameter (for which there is experimental information), while the individual 
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decays have an uncertainty from the decay constants (e.g. fs s ), and hence a reliance on 
theoretical techniques such as lattice QCD. Consequently, the double ratios of leptonic 
decays are an alternative and competitive probe of the parameter space of supersymmetric 
models. 
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